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ABSTRACT 


A  computer  program,  RADSPHERE,  has  been  developed 
for  calculating  the  axisymmetrlc  forced  vibration  response 
and  sound  radiation  of  a  thin,  elastic,  spherical  shell, 
submerged  in  an  infinite  fluid.  Four  cases  of  excitation 
are  provided:  (1)  a  concentrated  radial  load  acting  at  a 
point  on  the  shell  surface,  (2)  two  such  loads  acting 
simultaneously  at  antipodal  surface  points,  (3)  a  uniform 
load  distributed  over  a  sector  of  the  surface,  and  (4) 
two  sector  loads  acting  simultaneously  with  respect  to  an 
axis  of  the  sphere.  The  calculated  quantitites  are  the 
shell  surface  radial  velocity  and  the  acoustic  pressure 
at  both  the  shell  surface  and  in  the  far-field.  RADSPHERE 
results  are  generated  in  tabular  form  and  in  graphic  form, 
the  latter  through  use  of  a  companion  program  PLOTTER 
which  provides  an  interface  with  the  facilities  of  the 
software  package  DISSPLA. 
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This  work  was  sponsored  by  the  Foundation  Acoustic  Design  Program  at  the  David 
Taylor  Naval  Ship  Research  and  Development  Center  under  Task  Area  SI255SL00I,  Program 
Element  63569N,  and  Work  Unit  1211-601.  The  Naval  Sea  Systems  Command  program 
managers  were  Owen  Ritter  and  Richard  Chiu  (NAVSEA  Code  55Y3) . 


MATHEMATICAL  FORMULATION 

RADSPHERE  was  developed  in  response  to  the  need  to  calculate  analytic  solutions 

for  validating  a  finite  element  capability  called  NASHUA  ( NAS t ran  Helmholtz  Underwater 

Acoustics,  Everstine  et  al.S*  currently  under  development. 

The  formulations  which  the  program  evaluates  are  taken  either  directly  or  are 

derived  from  the  analysis  of  axisymmetric  forced  vibration  of  submerged  thin 

2 

spherical  shells  by  Junger  and  Feit  .  To  facilitate  referencing,  exact  mathematical 
expressions  taken  from  this  source  will  be  designated  by  the  notations  J  &  F  with  the 
original  equation  numbers,  in  addition  to  the  numbering  required  in  the  present 


report. 


A  convenient  starting  point  for  the  analysis  to  be  presented  here  is  the 
expression  for  forced  radial  velocity  response  of  the  shell  in  vacuo  (J&F,  Eq.  10.10) 


(cos6) 


where  6  -  the  colatitude  or  polar  angle  with  respect  to  an  axis  of  the  sphere 

Pn(co80)  ■  Legendre  polynomial  of  order  n  and  degree  1 

■  in-vacuo  modal  Impedance  of  the  shell 

f  •  coefficients  in  the  Legendre  series  expansion  for  an  applied  load 

acting  on  the  shell. 

The  in-vacuo  modal  impedance  is  given  by  (J&F,  Eq.  10.11), 


C02-(0^1))21  [fi2-(^2))2) 

_  .  n  n  n 

Z  -  -i-  p  c  - = - = - 

n  asp  n3  2.  ,  ,  . 

fi  -0(n  +n-l+v) 


where  i  *  /-T 


h  ■  shell  thickness 


a  -  shell  radius 


Pa  -  mass  density  of  the  shell  material 


v  “  Poisson' 8  ratio 


with  E  denoting  Young's  modulus 


0  ■  dimensionless  frequency  of  vibration 


•ft) 


ka  with  u>  ■  frequency  (rad/s) 
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c  ■  speed  of  sound  in  the  fluid 
k  ■  wave  number  ”  u/c 


0^^,  ”  dimensionless  frequencies  of  the  in-vacuo  shell  lower  and  upper 

n  n 

branch  modes  n,  respectively. 

The  total  loading  on  the  shell  is  expressed  as  a  combination  of  the  driving  force 
and  radiation  load  as 


p(e)  -  p(e)  +  p(e). 

total  applied  rad 


For  axisymmetric  modes  (J&F,  Eq.  10.1) 


p(6) 

rad 


-z 


z  W  P 
n  n  n 


(cos  6)  . 


(4) 


Here  z^  is  the  specific  acoustic  impedance  given  by 


h  (ka) 
,  n 

pfluidc  h' (ka) 
n 


“h*"  OfXuld 

h  (ka) 
n 


w 

n 


mass  density  of  fluid 

a  spherical  Hankel  function  of  first  kind  and  order  n  with 
argument  ka, 

the  modal  velocity  amplitude. 


The  analysis  now  requires  that  the  applied  load  p(6)  be  expanded  in  a 

applied 

Legendre  series  which  can  then  be  combined  with  the  series  in  Eq.  4,  to  yield  the 
coefficients  fR  of  the  expansion  of  the  total  load 


total 


n“o 


f  P  (cos  0) 
n  n 


( 


POINT  LOAD 

The  first  applied  loading  to  be  considered  here  is  illustrated  in  Fig.  1. 

In  this  figure  Fq  denotes  the  force  magnitude  and  t  is  the  time  variable.  The 

direction  of  the  force  indicated  is  positive  outward  toward  the  fluid.  The  force 

3 

is  represented  with  the  Dirac  delta  function 


F  (cos  0) 


F 

o 


6 (cos  6-1) 
2na2 


0 


Let  the  Legendre  expansion  of  this  force  be  represented  by 


F  (cos  6)  ■  /  f  P  (cos  0) 

i  »  m  m 

m 

where  f-m  denotes  the  expansion  coefficients. 
Then, 


F  (cos  0)  P  (cos  0) 
n 


(cos  0) 


P  (cos  0). 
n 


0  =  0° 

h 


( 


(I 


Fig.  1.  Concentrated  load  acting  on  a  submerged  spherical  shell. 


Integrating  Eq.  8  yields 


/F  (cos  9)  P  (cos  6)  d  (cos  9)  ■  f  ^  ^  f  P  (cos  9)  P  (cos 
n  J,  i  «  mm  n 


P  (cos  9)  P  (cos  9)  d 
m  n 


-  f 


2 

n  2n+l 


utilizing  the 


orthogonality  relations  given  by  Wylie 


4 


(COB  0)  P 

n 


(cos  9)  d  (cos  9) 


m^n 

m»n 


Returning  to  the  left  side  of  Eq.  9  there  results 


F  (cos  0)  P  (cos  9)  d  (cos  0) 
n 


6(cos  0  -  1)  P  (cos 
n 


p  a) 


2iTa 


2  n 


2na 


d  (cos  0) 


(cos  0) 


(9) 


(10) 


0)  d  (cos  0 


(ID 


4 

when  utilizing  the  fact  that  , 


f  »**_#■  f 


,-fc  ( 


*4,*  l,* 


f(x)  6(x-x  )  dx  -  f(x  ) 
o  o 


and  the  identity  PnC 1 )  *  1  for  all  n. 


From  Eqs.  9  and  11  one  obtains 


n  2n+l  _  2 

2ira 


from  which 


(2n+l)F 


n  4ira2  ’ 


Combining  the  coefficients  /  with  those  of  Eq.  4  yields  the  coefficients  fR  of  the 


total  load  for  this  case  (J&F,  Eq.  10.9) 


(2n+l)  F 


n  ,2 
4wa 


-  z  W  . 
n  n 


In  view  of  the  fact  that  in  Eq.  1  f  /Z  -  W  ,  one  has  by  substituting  for  f  from 

n  n  n  n 


Eq.  15 


(2n+l)F 


-  z  W 
2  n  n 


ni^nn7!nZEnRiSCVXRXi!CRni 


This  can  be  solved  for  (J&F,  Eq.  10.12),  obtaining 


(2n+l)  Fo 

4ira^(Z  +z  ) 
n  n 


(17) 


Substituting  Eq.  17  into  Eq.  1  gives  the  shell  radial  velocity  for  the  submerged 
case.  (J&F,  Eq.  10.13) 


w(8) 


(2n+l) 


Z  +z 
n  n 


P 

n 


(cos  6). 


(18) 
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Following  Junger  and  Felt's  analysis,  Eq.  18  is  then  multipled  by 
ip^^u^<jChn(kR)/h^(ka)  to  obtain  the  resulting  sound  pressure  for  field  points  at 
(R,  0)  where  R  is  the  radius  from  the  center  of  the  sphere.  If  R-a,  the  multiplying 
factor  becomes  z^,  Eq.  4,  and  one  has  the  expression 

F  (2n+l)  z 

P(a.e)  “  Lu~TTz —  Pn  (cos  6)*  (19) 

4na  n  n  n 

for  the  shell  surface  pressure. 

For  the  far-field  pressure,  hn(kR)  in  the  multiplying  factor  is  approximated  (J&F, 

Eq.  8.11c)  by 


It  is  readily  verified  that 


(20) 


7 


(21) 


Hence,  the  multiplying  factor  is  approximated  by 


P 


fluid 


(-i)n 


(22) 


and  the  far-field  formulation  corresponding  to  Eq.  19  is  then  (J&F  Eq.  10.14) 


p(R,0) 


F  p,,  .,c  e 
o  fluid 

4xa^kR 


ikR 


L(-i)P(2n+l) 

(Z+z)  h '  (lea)  *n 
n  n  n  n 


(cos  0) ,  kR  >> 


1. 


(23) 


Equations  (18,  19,  and  23  are  the  formulations  evaluated  by  RADSPHERE  for  the 
case  of  a  concentrated  harmonic  force  acting  at  the  apex  of  the  polar  axis. 

SECTOR  LOAD 

The  second  applied  loading  of  interest  is  illustrated  in  Fig.  2.  This  load 
is  defined  by 


{F  ,  o  <  6  <  a  <  tt 
° 

o,  otherwise. 


(24) 


For  this  case  the  left  side  of  Eq.  9  is 


/ 


F  (cos  0)  P  (cos  0)  d  (cos  0)  ■  F 
n  o 


/  Pn 


(cos  0)  d  (cos  0).  (25) 


cos  a 
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IT 

K 

[< 


0-0° 


b/. 


UNIFORM  LOAD: 

F  -  K  •-*"» 


Fig.  2.  Uniform  load  acting  over  a  polar  sector  of  a  submerged  spherical  shell. 


Utilizing  the  identity*. 


pn+l(x)  “  Pn-l(x)  ’  (2n+1)  Pn(x)’ 


Eq.  25  yields 


Fo  J  Pn(cos  6)  d(co8  Q)  ~  Z^lJ  {P 
•cob  a  »cosa 


' , , (cos  9)-P'  .(cos  0)]  d(cos  0) 
n+i  n— i 


2S+I  (p„+l(cos  9)-p„-l(cos  6)Jcos- 


2STT  u-‘-[pn+i(cos  “)-p„-i<co»  “>» 


25+T  [Pn-l(co8  o)-Pn+l(co8  0)1 


with  the  Identity  for  Pn(l)  given  on  page  6. 

Equating  this  result  to  the  right  side  of  Eq.  9  and  solving  for  f-  leads  to 


(27) 


fnmT  lPn-l  (cos  o)  "  Pn+1  (c08  “)J  * 


(28) 


Analogous  to  the  steps  taken  in  Eqs.  15-17,  respectively,  one  obtains  the  total  load 


f  “  [P  ,  (cos  a)-P  (cos  a) ]  -  z  W 
n  i  n— i  n+i  n  n 


(29) 


then. 


~2  tpn_i(cos  a)“pn+l (cos  a)J 


z  W 
n  n 


n 


and,  solving  for  W  , 
n 


(30) 


W 


~2  fpn_i (cos  a)-pn+i (c°s  a)] 

(Z  +  z  ) 
n  n 


(31) 


P 


Substituting  for  (fn/Zn)  in  Eq. 


» 

i 

1  gives  the  shell  radial  velocity 


(32) 


w(e) 


p  v**£p 

-fL— 


n*o 


1  (cos  a)  "Pn+1  (cos  a) ] 

<z„  +  *  > 

n  n 


p  (cos  e) 
n 


from  which  is  obtained,  utilizing  steps  analogous  to  those  used  for  Eqs.  19  and  23 
respectively,  the  shell  surface  pressure 


p(a 


F  IP  . 

.»  --fL  — 


(cos  a)  -  Pn+1  (cos  a)]  z r 


n“o 


(Z  +  z  ) 
n  n 


P  (cos 
n 


0) 


(33) 


and  the  far-field  pressure 


p(R.e)  ■ 


F  pri  ,  ,ce 
o  fluid 

2kR 


lkR 


n*o 


°>  -  Vi 
(2»  +  *„>  h;<i“) 


(cos  a) ] 


P  (cos 
n 


6).  (34) 


Equations  32-34  are  the  formulations  evaluated  RADSPHERE  for  the  case  of  an 
harmonic  driving  force  distributed  uniformly  over  a  sector  of  the  shell  surface 
centered  on  an  apex  of  the  polar  axis. 

The  shell  response  and/or  field  pressure  for  two  additional  load  conditions  is 
readily  obtained  by  superposing  solutions  for  the  loading  conditions  just  discussed. 
This  process  is  pictured  in  Fig.  3.  If  a  complex-valued  solution  corresponding  to 
load  condition  I  is  designated  Solj  (e)  and  for  load  condition  II  by  Sol^  (0), 
then 


SolII  (0i)  "  Soll  ^18°-ei) 


(35) 


and 


11 


k 


Fig.  3.  Symmetric  loads  synthesized  from  asymmetric  load  components 
(denoted  as  load  conditions  I  or  II) . 


Solsym  (0i)  "  Soli  (0i>  +  Soli  O8O-0i)  -  Solj  (0^  +  Soln  (0^. 


The  basic  calculations  made  by  RADSPHERE  are  for  load  conditions  I,  which  are 


designated  as  PROBLEM  1  types  in  the  program  data  input.  The  two  symmetric  conditions 


are  designated  as  PROBLEM  2  types. 


NOTES  ON  THE  COMPUTATION 


I.  Calculation  of  the  specific  acoustic  impedance  (Eq.  4)  requires  spherical 


Hankel  functions  and  their  derivatives  for  particular  arguments.  In  standard 


notation,  the  spherical  Hankel  function  h  of  order  n  for  argument  x  is  represented 


K 


t- * 


$ 


$ 

& 


R 


i 


6: 


w: 


$ 

$ 


5! 


« 

5S 


(37) 


h  (x)  -  j  (*)  +  iy  (x) 
n  n  n 


where  1  (x)  and  y  (x)  are  Che  spherical  Bessel  functions  of  the  first  and  second 
n  n 


kinds,  respectively,  and  of  order  n.  To  obtain  these  latter  two  functions,  IMSl 
subroutines  MMBSJR  and  MMBSYN  are  first  used  to  calculate  ^n+j^/2^x^  an<*  Yn+1/2^X^ 


which  are  the  Bessel  functions  of  the  first  and  second  kinds,  respectively,  and  of 


fractional  order.  The  spherical  Hankel  functions  are  then  calculated  (Abramowltz 


and  Stegun  )  by  the  rule 


hn(x)  -  /tt/2x  lJn+1/2M  +  1Yn+l/2^X^'  •• 


(38) 


Derivatives  of  these  functions  are  obtained  by  the  rule 


h^(x)  “  2n+T  *nhn-l(x)  "  (n+l)hn+I (x) J ,  n-0,1,2,  .. 


(39) 


2.  Evaluation  of  the  Legendre  polynomials  of  order  n  and  degree  v«l,  for 

7 


argument  x,  are  obtained  by  the  recursion 


PQ(x)  -  1, 


Pn+1 (x)  •  x  pn(x) 


*  Pn-l(x> 


(40) 


3.  The  RADSPHERE  program  requires  specification  of  the  maximum  number  of  terms 
that  will  be  included  in  any  series  being  summed  (surface  pressure,  surface  velocity, 
field  pressure)  for  a  particular  application.  The  purpose  of  this  is  to  bound  the 
computing  effort  that  will  be  expended  trying  to  obtain  convergence  in  the  event  of 
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very  raall  magnitudes  (i.e.,  nodes  in  8  profiles  for  velocity  or  pressure).  Belov 
this  maximum  number*  the  series  summation  is  terminated  when  the  condition 


Sum 


iH til 


-  Sum 


ls“J 


<  e 


(41) 


is  met*  with  c  being  a  specified  error  bound.  The  index  n  references  only  nonzero 
terms  in  the  series  being  summed. 

4.  The  RADSPHERE  program  treats  the  frequency  range  of  interest  (in  terms  of 
ka)  as  divided  into  subintervals  i,  not  necessarily  adjacent,  each  having  uniform 
spacing  (Aka)  within.  This  treatment  facilitates  computing  for  a  set  of  nonunlformly 
spaced  ka's  (since  an  "interval"  can  consist  of  a  single  ka)  as  well  as  computing 
at  varying  densities  of  ka  over  the  range  of  interest.  The  latter  capability  is 
generally  useful  in  the  case  of  searching  for,  or  refining  resonant  frequencies  of 
the  submerged  shell. 


DATA  INPUT  TO  RADSPHERE 

g 

Data  for  RADSPHERE  is  input  with  list-directed  (unformatted)  records.  Some 
of  the  data  control  option  selections  and/or  branching  within  the  program.  These 
parameters  are  designated  by  their  FORTRAN  variable  names.  The  data  are  presented 
in  Table  I  in  the  order  required  by  the  program.  In  this  table  and  all  subsequent 
tables  throughout  the  report*  each  line  in  the  column  headed  "Data*'  represents  one 


card  record  of  data. 


Title 


NFLAGF 


pfluid 


NDIAGN 


NPROB 


Table  1.  RADSPHERE  data  input. 


Description 


50  alphanuaeric 
characters  in  card 
columns  1-50 


Load  selection 


Hagnitude  of  load 


Polar  angle  (deg)  of 
sector  load  (Fig.  2) 


Shell  radius 


Shell  thickness 


Shell  mass  density 


Poisson's  ratio 


Young's  modulus 


Structural  loss  factor 


Fluid  mass  density 


Speed  of  sound  in  the 
fluid 


Controls  printing  of 
intermediate  quantities 
for  check  purposes 


Selector  for  printing 
calculated  results 


Application 


When  NFLAGF-2 


Value 


1  (point  drive  at 
6-0"  (Fig.  D) 


2  (sector  drive 
(Fig.  2)) 


0,  or  a  specified 
value 


1  (Do  not  print) 

2  (Print) 


1  (Problem  type  1 
(Fig.  3)  results 
only) 


2  (Problem  type  2 
(Fig.  3)  results 
only) 


3  (Both  problem 
types  1  and  2, 
results) 


m 


,  v.V/v 


Table  1.  (Continued) 


Data 

Description 

Application 

Value 

NFLAG 

Selector  for  quantities  to  be 
calculated 

■  1  (shell  surface 
pressure) 

-  2  (shell  surface  and 
field  pressure) 

■  3  (field  pressure) 

■  4  (shell  surface 

pressure  and  radial 
velocity) 

NRADII 

Number  of  radii  at  which  field 
calculations  are  to  be  made, 
maximum  -  10 

When 

NFLAG  - 
2  or  3 

■  5  (shell  surface 
radial  velocity) 

l' 

2 

• 

• 

Distinct  radii  from  center 
of  shell  to  field  points. 

When 

NFLAG  - 
2  or  3 

• 

Sradii 

NTERMS 

Bound  on  no.  of  terms  In  any 
series  to  be  computed  (see 
pp.  13-14),  maximum  ■  200 

All 

e 

Error  criterion  for 
convergence  of  series 
summations. 

All 

61 

Initial  colatitude  (deg)  for 
shell  surface  and/or  field 
points  at  which  calculations 
are  to  be  made  (see  Fig.  3). 

All 

02 

Final  colatitude  (deg) 

All 

A0 

Increment  in  colatitude  (deg) 

All 

[  4,  1  +  *]  *  91 

NINT 

Number  of  intervals  in  ka 
(see  page  14) 

All 

Data 


Table  1.  (Continued) 


Data 

Description 

Application 

Value 

(kajJj,  (ka2>1,  (Aka) j 

^^*2^2*  (Alcfi) 2 

•  •  • 

•  •  • 

Initial  ka,  final 
ka,  Aka  for  the 
i^  interval  in  ka. 

All 

•  •  • 

^l^NINT*  (k^NINT*  (Aka) jjijjT 

With  regard  to  the  data  pertaining  to  the  shell  geometry  and  the  shell  material 
and  fluid  properties,  it  is  noted  that  any  system  of  units  may  be  used  so  long  as 
consistency  is  maintained. 


The 

data  flag  NDIAGN  is 

used  to  obtain  printout  of  the 

following  quantities  for 

diagnostic  purposes: 

k®]1 

(see  p.  3) 

b. 

Pn  (cos  a) 

(see  Eq.  28) 

c. 

Jn+l/2  (ka)’  Yn+l/2 

(ka) 

(see  p.  13) 

d. 

hn  (ka) 

(see  Eq.  38) 

e. 

h'  (ka) 
n 

(see  Eq.  39) 

f . 

z 

n 

(see  p.  3) 

g- 

Z 

n 

(see  p.  2) 

h. 

z 

n 

(see  Eq.  19) 

(Z+z) 
n  n 

i. 

1 

(see  Eq.  23) 

(Z+O  h'  (ka) 
n  n  n 

J. 

1 

(see  Eq.  18) 

i 

(Z  +8  ) 
n  n 

k. 

P  (cos  0) 
n 

(see  Eq.  1) 
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The  NDIAGN  is  also  used  Co  echo  the  data  written  to  plot  files. 

To  prevent  the  "wallpaper"  effect  when  using  NDIAGN,  it  is  advisable  to  trim 
diagnostic  runs  to  the  minimum  number  of  ka's  (ideally,  -  1)  and  as  low  a  value  of 
"NTERMS"  as  possible. 

The  loss  factor  n  is  introduced  into  the  calculations  via  the  familiar  mechanism^ 
of  scaling  E  by  (1-in).  The  sign  of  in  is  consistent  with  the  exponential  form  for 
the  time  function  e  *u>t  assumed  in  the  analysis  by  Junger  and  Feit.^ 

The  RADSPHERE  program  is  dimensioned  for  a  maximum  of  91  values  of  colatitude 
for  computations  of  response  profiles,  and  a  maximum  of  2001  values  of  ka  per 
interval  of  ka.  The  number  of  ka  intervals  is  arbitrary. 

For  calculations  involving  only  a  single  colatitude,  NPROB  must  be  assigned  the 
value  of  one.  This  restriction  results  from  the  fact  that  the  single  colatitude 
case  was  designed  exclusively  to  facilitate  searching  for  resonant  ka.  While  either 
single  or  double  pole  excitation  could  be  used  for  this  case  (see  Fig.  3),  only  the 
single  excitation  response  is  efficient  to  calculate  because  only  a  single  solution 
is  required  (see  Eq.  36). 

DATA  INPUT  TO  PLOTTER 

Certain  subsets  of  the  results  computed  by  RADSPHERE  can  be  obtained  in  graph 
form  via  the  companion  program,  PLOTTER,  which  utilizes  the  software  package 
DISSPLA^.  In  particular,  one  can  obtain  polar  plots  of  shell  surface  pressure  and 
radial  velocity,  as  well  as  far-field  pressure,  for  each  ka,  or  log-linear  plots  of 
the  same  quantities  at  a  single  colatitude  over  a  range  of  ka  values. 

Data  to  be  plotted  are  written  by  RADSPHERE  on  a  permanent  file  device  TAPE12. 
Following  a  RADSPHERE  run  this  file  must  be  cataloged  for  subsequent  use  by 
PLOTTER  when  plots  are  desired.  Before  data  intended  for  polar  plots  are  written  on 


TAPE12,  a  subroutine,  SORT,  recasts  the  data  from  the  system  used  in  computing 
(see  Fig.  3)  to  the  coordinate  system  used  for  plotting.  Fig.  4.  In  the  plotting 
coordinate  system,  a  result  computed  at  colatitude  3-0°  appears  at  8  «  90°;  a 
result  at  colatitude  6  *  90°  appears  at  0  *  0°;  a  result  at  colatitude  6  »  180° 
appears  at  8  •  270°;  etc.  Subroutine  SORT  also  augments  (via  axial  symmetry)  the 
original  data  computed  in  the  colatitude  range  0°  <  9  <_  180°  to  the  full  range  in 
plotting  coordinates,  0°  <  0  £  360°. 

For  the  variables  shell  surface  pressure  and  field  pressure,  the  quantity 
plotted  is  always  the  absolute  value.  For  log-linear  plots  of  shell  surface  velocity 
versus  ka,  the  plotted  quantity  is  also  the  absolute  value.  However,  for  polar  plots 
of  velocity  profile,  it  seemed  advisable  to  attempt  to  Incorporate  the  effect  of 
phase.  An  approach  is  to  plot 

¥ 

w(0)  -  sin  y  •  | w(0)  |  -  J  w(0),  (42) 

where  y  ■  tan 

L^(0)  J 


270° 


Fig.  4.  Coordinate  system  for  plotting. 


with*/,  ^designating,  respectively,  imaginary  and  real  parts  of  w(0).  The  choice 
of  sin  y  was  made  to  obtain  relative  maxima  of  response  at  0  *  90°  and  9  ■  270°, 
as  induced  by  the  excitation  force(s). 

Much  of  the  input  data  required  by  PLOTTER  (aside  from  the  plot  data)  is 
supplied  on  the  plot  file  from  RADSPHERE.  The  remaining  input  data  required  is 
largely  that  necessary  to  provide  the  user  considerable  flexibility  in  shaping  the 
plots. 

The  following  set  of  data  (Table  2)  is  provided  for  each  interval  of  ka 
computed. 


Table  2.  Plotter  data  input. 


Application 


NCYCLE,  YMIN 


XSTEP 


NXTICK.NYTICK 


Log-linear 
plots . 


Table  2.  (Continued) 


Description 


Shell  Surface  Velocity 
Problem  type  1 
Problem  type  2 

Number  of  cycles  and  minimum  value  for  the  y 
axis. 


Step  size  in  the  x-direction.  Note:  x  ^  + 

10*XSTEP  -  x 


Parameters  controlling  the  placement  of  grid 
lines  in  the  x-  and  y-directions ,  respectively. 
For  the  purpose  here,  the  particular  form  of 
these  parameters  is  NXTICK  *  -m,  NYTICK  *  -n, 
which  places  grid  lines  at  every  mc^  and  nc^ 
subdivision  on  the  respective  axes. 

Note:  This  group  of  three  cards  is  provided  in 
each  ka  interval  for  each  computed  variable  in 
the  same  order  as  listed  above  for  polar  plots. 
Only  type  1  problem  data  will  be  present  as 
previously  stated. 


All  data  for  PLOTTER  is  specified  in  list-directed  form  . 

The  output  from  PLOTTER  consists  of  a  single  page  of  printed  output  which  echos 
identifying  data  from  the  RADSPHERE  data  file,  and  PLF1LE,  written  by  DISSPLA.  The 
latter  file  is  cataloged  for  postprocessing  either  to  a  Tektronix  screen,  inter¬ 
actively,  or  to  a  tape  for  off-line  Calcomp  plots. 


CALCULATIONS 

The  final  section  of  this  report  presents  several  calculations  to  illustrate 
some  of  the  features  of  RADSPHERE  that  have  been  discussed  and  to  provide  examples 
of  the  various  forms  of  edited  output. 


“ y  v*>y«w 
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UNIFORMLY  PULSATING  SPHERICAL  SHELL 


For  this  problem  (see  Fig.  5),  the  submerged  shell  is  driven  by  a  uniformly 
distributed  pressure  of  magnitude  Fq  acting  with  harmonic  time  variance  over  the 
interior  shell  surface.  This  is  an  ideal  problem  from  the  standpoint  of  validating 
RADSPHERE,  because  closed  form  solutions  can  be  derived^  for  the  shell  surface 
pressure  and  radial  velocity  and  the  field  pressure.  Figures  2  and  3  show  that  this 
problem  can  be  cast  in  the  form  of  superposed  loadings  in  which  I  and  II  are  sector 
loads  with  polar  angle,  a  *  90°.  In  RADSPHERE  terminology,  this  is  a  Problem  2 
type  example.  The  calculation  is  performed  in  two  runs:  a.  to  obtain  the  surface 
pressure  and  radial  velocity,  and  b.  to  obtain  field  pressure.  The  data  input  is 
shown  in  Table  3.'  Data  for  the  shell  are  in  MKS  units  for  steel. 


Fig.  5.  Submerged  uniformly  pulsating  spherical  shell. 
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Table  3.  RADSPHERE  input  data  for  problem  of  uniformly 
pulsating  spherical  shell. 


Problem  Part 

Data  Item 

Value 

Part  a,  to  calculate 
surface  pressure  and 
radial  velocity. 

Title 

NFLAGF 

Uniformly  Pulsating  Sphere 

2 

F 

o 

1. 

a 

90. 

a 

5. 

h 

0.15 

ps 

7669. 

V 

0.3 

E 

2.07E11 

n 

0. 

p fluid 

1000. 

c 

1524. 

NDIAGN 

1 

NPROB 

2 

NFLAG 

4 

NTERMS 

30 

e 

0.0001 

01 

0. 

e2 

180. 

A0 

10. 

NINT 

1 

(kaj )1 

,  (ka2)1>  (Aka)L 

5.,  5.,  1. 

Part  b,  to  calculate 
field  pressure 

Title 

• 

same  as  for  a. 

NPROB  1 

NFLAG 

3 

NRADII 

1 

R1 

100. 
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Table  3.  (Continued) 


Problem  Part 

Data  Item 

Value 

NTERMS 

• 

• 

• 

(ka1)1,  (ka2)1,  (Aka)1 

same  as  for  a. 

The  respective  printouts  from  RADSPHERE  for  these  sets  of  data  are  shown  In 
Fig.  6  (Part  a)  and  Fig.  7  (Part  b) .  Referring  to  the  tabular  output  in  these 
figures,  the  columns  headed  ERROR  1  and  ERROR  2  give  the  error  achieved  at 
convergence  for  components  I  and  II,  respectively  (see  Fig.  3)  of  the  solution. 
The  columns  headed  FLAG  1  and  FLAG  2  give  the  number  of  series  terms  required 
to  achieve  the  respective  levels  of  convergence  and  thus  are  indicators  of  the 
degree  of  difficulty  in  obtaining  the  solution  at  each  colatitude. 


PRORLEN  TITLE  -  UNI  FOR  Hi T  PULSATING  SPHERE 


DATA  FOR  SHELL 


RADIUS  >9*00 

THICKNESS  •  .19 

HASS  DENSITY  -T.MD0C403 

POISSON  RATIO  •  .30 

TOWNS  S  NOOULUS-2.0TOOE411 
LOSS  FACTOR  >0*00 


DATA  FOR  EXCITATION 
SECTOR  LOAD 

POLAR  ANDLE  •  90.00  OESREES 
HASNITUOE  •  1.00 


DATA  FOR  FLUIO 

HASS  OENSITY  •1000.00 
SPEEO  OF  SOUNO  >1924.00 


PROSLEH  TYPES  TO  SE  CALCULATED 


2 


QUANTITIES  TO  RE  COHPUTEO 


SHELL  SURFACE 
VELOCITY 


SHELL  SURFACE 
PRESSURE 


X 


X 


DIAGNOSTICS 


00  NOT  PRINT 


CALCULATIONS  TO  RE  NAOE  FOR  - 


INITIAL  COLATITUOE  -  0.00 
FINAL  COLATITUOE  *1Q0.00 
OELTA  COLATITUOE  •  10.00 


NAX.  NO.  OF  TERNS  FOR  A  SERIES  ■  SO 


C0NVER6ENCE  CRITERION  •  1.0000E-04 


KA  INTERVAL  NUHSER  IS  1 

INITIAL  KA  ■  5 .00001 ♦00 
PINAL  KA  •  5, 00001 ♦ 00 
DELTA  KA  •  1.00001400 


FAR  -  FIELD 
PRESSURE 


Fig.  6.  RADSPHERE  output  for  uniformly  pulsating  spherical  shell  (Part  a) 


MONISH  TITLE  •  UNIFORHLT  FULSATIN6  SPHERE 

MONLSR  2.  -  OUTWARD  OISTRIRUTED  FORCE  AT  ROTH  NORTH  ANO  SOUTH  ROLES* 
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FR0RLEN  TITlt  -  UNIFORNLY  PULS AT INC  SFHERE 


MTI  roc  SHELL 

RAOIUS  *9.00 

THICKNESS  -  *19 

NASS  OINSITY  -T.66001P03 

POISSON  RATIO  •  .30 

VOUNS  S  NODULU$«2.0700E*11 
LOSS  FACTOR  -0*00 


rata  for  excitation 


SECTOR  LOAO 

FOiAR  ANCLE  •  CO.OO 

nacnituoe  •  1*00 


OECRECS 


RATA  FOR  FUIO 

NASS  DENSITY  *1000.00 
SPEED  OF  SOONO  *1924.00 

RATA  FOR  FIELO 

NO.  OF  RAOII  •  1 

RAO IUS(  1)  *100.00 
FRORLIN  TYFES  TO  RE  CALCULATED 

2 


QUANTITIES  TO  DC  CONFUTED 

SNELL  SURFACE 
VELOCITY 


SHELL  SURFACE 
FRESSURE 


FAR  -  FIELD 
FRESSURE 

X 


R1ACN0STICS 

00  NOT  FEINT 

CALCULATIONS  TO  RE  NAOE  FOR  - 

INITIAL  COLATITUOC  •  0.00 

FINAL  COLATITUOi  *110.00 
REITA  COLATITUOi  •  10.00 


Pig.  7.  RADSPHERE  output  for  uniformly  pulsating  spherical  shell  (Part  b) 
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MAX.  NO*  OF  TERNS  FOR  A  SERIES  ■  30 


CONVERGENCE  CRITERION  ■  l.OOOOE-OA 


INTERVAL  NUNRCR  IS  1 

INITIAL  M  •  9*OOOOE«0O 
FINAL  FA  •  S.COOOEtOO 
PELTA  FA  -  l.OOOOEROO 


Fig.  7.  (Continued) 


30 


PROBIEN  T1TIC  -  UN  I  FOR  NIT  POLSATINB  SPHERE 

P ROUEN  2.  -  OUTRAN?  DISTRIBUTED  FORCE  AT  DOTH  NORTH  AND  SOUTH  ROLES* 
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Table  4  shows  the  numerical  solution  to  be  In  excellent  agreement  with 
analytic  results^. 


Table  4.  Comparison  of  RADSPHERE  and  analytic  results  for  the 
uniformly  pulsating  spherical  shell. 


RADSPHERE 

Analytic 

Shell  Surface 

6.5800E-7  -  il . 2674E-7 

6.58E-7  +  il . 27E-7 

Radial  Velocity 

Shell  Surface 

9.2708E-1  -  13. 7857E-1 

9.27E-1  +  13.79E-1 

Pressure 

Field  Pressure 

4.6780E-2  +  il . 7851E-2 

4.68E-2  -  il . 79E-2 

The  difference  in  sign  of  the  imaginary  parts  reflects  the  fact  that  the  time 
dependence  for  RADSPHERE' s  analysis  is  e  iwt,  whereas  that  for  the  analytic 
solution  is  eiait. 

POINT-DRIVEN  SPHERICAL  SHELL 

The  second  calculation  is  for  a  spherical  shell  whose  material  and  geometric 

3 

properties  and  excitation  force  are  taken  from  Hayek's  investigation.  The 
physical  problem  represented  here  is  that  shown  previously  in  Fig.  1.  The  shell 
surface  pressure  and  radial  velocity  are  computed,  and  the  results  graphed  in  polar 
form  with  PLOTTER.  The  Input  data  for  RADSPHERE  is  given  in  Table  5. 

3 

The  value  of  ka  is  derived  from  an  expression  for  the  frequency  parameter  , 


2  2  2 
fi  *  (1/E)  p  o>  a 
s 

with  £)  given  the  value  of  0.85.  Data  for  the  shell  are  in  in-lb-s  units  for 

steel. 


Z+Z+Z&'CZ’t&i 


Table  5.  RADSPHERE  Input  data  for  problem  of  point-driven  spherical  shell. 


Data  Item 

Value 

Title 

Point-Driven  Sphere 

NFLAGF 

1 

F 

1. 

o 

a 

5. 

h 

0.15 

P 

0. 7347E-3 

s 

V 

0.3 

E 

30.  E6 

n 

0. 

p,.  .  . 

9.580488E-5 

fluid 

c 

60000. 

NDIAGN 

1 

NPROB 

1 

NFLAG 

4 

NTERMS 

90 

E 

0.0001 

91 

0. 

02 

180. 

A6 

2. 

NINT 

1 

(ka^) t ,  (ka2)1>  (Aka>1 

2.86268316,  2.86268316,  1. 

RADSPHERE  output  is  shown  in  Fig.  8.  Obviously,  from  the  number  of  terms 
required  (column  "FLAG"),  obtaining  convergence  of  velocity  at  the  poles  is 
considerably  more  difficult  than  for  the  corresponsing  pressures. 


DATA 


ECHO 


RROJLFp  TITLE  -  ROINT-CAIVfN  SPHERE 


DATA  FO»  SMELL 

RADII'S  •  S.OO 

THICKNESS  -  .15 

NASS  DENSITY  •7.3470E-05 

POISSON  RATIO  •  .30 

YOUNG  S  N0DIIUS-7.0000EAC7 
LOSS  FACTOR  -e.oo 


DATA  FOR  EXCITATION 

POINT  LOAD 

“ACNITUOE  ■  l.OC 


DATA  FOR  FLUTO 

NASS  DENSITY  •Q.590ABBE-CS 
SPEFP  ne  SOUND  -6 .Ot COOOE *05 


•POPLEP  TYPES  TO  BE  CALCULATED 


1 


QUANTITIES  TO  Be  CONFUTED 

SHELL  SURFACE  SHELL  SURFACE  FAP  -  FIELD 

VELOCITY  PRESSURE  PRESSURE 

X  X 


01  AGNOSTICS 


DO  NOT  PRINT 


CALCULATIONS  TO  BE  NAOfc  EOS  - 


INITIAL  fOlATITUDE  •  0.00 
FINAL  COL  ATITUOF  *190.00 
DELTA  COLATITUOt  •  ?,00 


NAX.  NO.  OF  TF*NS  FOR  A  SERIES  •  00 


CONVERGENCE  CRITERION  •  l.cooOE-OS 


KA  INTERVAL  NUNREP  IS  1 

INITIAL  KA  •  ?.eP?7F»00 
FINAL  k*  ■  ?.«6?7E400 
DELTA  K  A  *  I.CflA-nE+l'l 


Fig.  8.  RADSPHERE  output  for  point-driven  spherical  shell. 


34 


uviNinininuvuvnivv 


mnumn 


.  NNICK  TITii  -  R01HT-0RIV6N  SNtfll 
PR00L6H  1.  -  OUTMARO  CONCENTRATED  FOKI  AT  MOUTH  POLE. 
Il<  2.0* 

FREOU6NCY  -946T.lt  HI 


C0LATITUDEI0E6R66S) 

PRESSURE 

AOSOIUTE  VALUE 

PHASEID66I 

ERROO 

FLAO 

0.00 

-3. 78*2 £-01 

6.91S2E-06 

I.T092E-01 

179.09 

4.998-09 

91 

I. 00 

-3, JT46E-01 

6.99648-04 

3.97448-01 

179.09 

9.438-09 

44 

4.00 

-3.20668-01 

7.24208-04 

3.20668-01 

179.87 

0.438-09 

12 

6.00 

— 2. 6707E-01 

7.44208-04 

2.47078-01 

179.04 

0.908-09 

23 

0.00 

“1.96001—01 

0.17748-04 

1.99008-01 

179.77 

9.00E-09 

37 

10.00 

-1.2004E-01 

0.02648-04 

1.20048-01 

179.61 

1. 44E-09 

32 

It  .00 

-9.3100F-02 

9.96148-04 

1.9188E-62 

179.01 

4.19E— 09 

27 

14.00 

1.344*1-02 

1.0392E-03 

1.34008-02 

4.40 

4.178-09 

36 

16.00 

T.3012E-O! 

1.11638-03 

7.9020E-02 

.04 

9.078-09 

42 

16.00 

1.2T34E-01 

1.19608-03 

1.27608-01 

.94 

1.298-09 

20 

to. 00 

1.6703E— 01 

1.27078-03 

1.67048-01 

.43 

2.418-09 

34 

tt.00 

1.4342E-01 

1.3369E-03 

1.9I42E-01 

.40 

1.018-09 

23 

14.00 

2.0610E-01 

1.39198-03 

2.06108-01 

.19 

4.076-09 

29 

26.00 

2.0304E— 01 

1.43198-03 

2.03048-01 

.40 

4.096-09 

99 

ta.oo 

l.OOOOE-Ol 

1.49448-03 

1.00818-01 

.44 

7.798-04 

11 

30.00 

1.61208-01 

1.49838-03 

1.6120E-01 

.92 

7.246-09 

IT 

32.00 

1 .232*8-01 

1.44108-03 

1.29306-01 

.66 

4.206-09 

14 

14.00 

0.464*8-02 

1.40398-03 

0.46616-02 

.99 

1.098-09 

31 

36.00 

3. 04918-02 

1.34448-03 

3.0914E-02 

2.00 

7.478-09 

39 

30.00 

-0.34038-03 

1.26368-03 

0.64128-03 

171.99 

4.696-09 

42 

40.00 

-9.36198-02 

1.16228-03 

9.36318-02 

170.76 

7.676-09 

91 

42.00 

-9.39338-02 

1.04198-03 

9.39308-02 

179.14 

6.908-09 

94 

44.00 

-1.27328-01 

9.03108-04 

1.27328-01 

179.99 

2.006-09 

24 

46.00 

-1.92148-01 

7.40978-04 

1.92148-01 

179.72 

4.928-06 

29 

Fig.  8.  (Continued) 
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40.1./ 

-1.0  70JE-01 

5,0131*— 04 

1.67030-01 

170.60 

7.710-05 

22 

jt.ao 

-i.»«  oof-n 

5.07560-05 

1.70000-01 

170.66 

6.200-05 

14 

it.ij 

-1.63360-11 

7.15670-05 

1.65550-01 

170.03 

5.120-05 

31 

-1.40670-91 

7.07750-03 

1.50670-01 

170.00 

4.710-06 

23 

56.44 

-1.75506-91 

-1.773*0-05 

1.75500-01 

161. Ct 

4.310-05 

10 

50.00 

—0.57701—97 

-3.77|60-v5 

0.37710-02 

160.73 

4.100-05 

14 

60.60 

—f .67610— 3? 

-3.77316-05 

5.67650-07 

160.51 

0. 067-05 

36 

62.49 

-l.T6*S£-<.'2 

-7.75136—05 

1.76600-47 

162.57 

2.700-05 

23 

05.00 

7.15000-97 

-0.60050—05 

7.13770-07 

337.30 

7.040-06 

42 

66 . 3C 

5. 00050-07 

-1.15076-vl 

5.60060-07 

336.66 

4.230-07 

36 

60.00 

0.05657-07 

-1.35360— Cl 

0.65750-47 

330.15 

5.460-03 

21 

70.00 

1. If 070-01 

-1.37070-03 

1.15030-01 

350,25 

2.630-06 

23 

77. CO 

1.53050— Cl 

-1.60016-03 

1.33656-01 

350.76 

3. 530-03 

30 

75.09 

1. 57651—01 

-1.65150-03 

1.52660-01 

350.76 

4.110-05 

20 

76.  V 

1 .5100  c— 01 

-7.00506-03 

1.51070-01 

350.10 

3.260-05 

26 

70.00 

1 .37737-01 

-7.15760-03 

1.32250-01 

350.67 

1. 420-wO 

23 

00.00 

1  .13730-91 

-7. 76176-03 

1.13770-01 

356.65 

3.760-03 

27 

6?.*. 

6.70J7C- |J 

-7.50620-03 

8.70500-02 

356.43 

5.060-05 

31 

05.01 

3.«503>-C7 

-?.S2C7t-C3 

1 5.65*«t-0? 

357.45 

4.200-03 

36 

66.'-< 

0. 176  5? — 9? 

-7. 6250* -03 

2.15270-02 

352.06 

2.630-06 

23 

00.00 

-l.S71«f-<.? 

-7.71666-51 

1.55500-02 

100.11 

5.060-03 

43 

i 

-«.<  0000-97 

-7.60166-03 

5.00150-0? 

163.15 

0.760-03 

40 

07.  v. 

-6.|503f-.'7 

-6.67366-03 

6.35140-02 

161.07 

0.000-0! 

33 

05.  at. 

-l.iCS-f-.l 

-7 .035*0— C  3 

1.10630-01 

161.52 

3.510-07 

21 

06.00 

-1. 71111-01 

-7.e*3*'-03 

1.11341-91 

161.30 

1. 140-03 

10 

oo.n.- 

-l.5333t-'l 

-3  «i'27ic-t'3 

1.43560-01 

181.21 

3.650-05 

33 

Uj." 

-1.56730-1 1 

-3.1/5016-^3 

1.46756-91 

161.10 

4.470-05 

27 

1C7.0C 

-1.5C«»E-C1 

-3,96517—03 

1.40616-41 

181.25 

1.330-06 

23 

195. 99 

-l.2534fc-<  l 

-3, .71 5f-i  3 

1.75300-61 

161.43 

3.710-03 

76 

lC6.*/'< 

-l.tlOOf-'  i 

-3. .6556-93 

1.02040-91 

161.72 

3.730-03 

20 

1C5.0? 

-7.1»7?o-07 

-l.e5O*0-6J 

7.10371-0? 

182.43 

3.100-03 

33 

11  J.  10 

-3.6  7510— C7 

-5.92756-03 

3.66660-07 

164.71 

8.800-96 

73 

117.00 

1.0695--.I 

-7,006ut-</3 

1,16760-43 

280.07 

6.650-05 

36 

« 

Fig.  8.  (Continued)  - 
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u*.  pe 

1.1*616-07 

-2.9*7*6-01 

1.99716-02 

393. Tl 

6.21E-07 

u».oo 

?•?  9306— 12 

-2.91136-03 

7. 39166-92 

137.10 

7.516-05 

111. 00 

1.06066-01 

-2,66*16— ol 

1.01116-01 

331.46 

4.37E-06 

Mi.io 

1 .1*0*6—01 

-2.*ll?«-03 

1.3*076-01 

356.60 

9.576-03 

Ift.w 

l.«162*-0l 

-2.7I566-C3 

1.31156-01 

336.96 

6.866-03 

It*. 90 

1.60196-01 

-2.69676-03 

1.60216-01 

339.03 

3.336-05 

m.*t 

1.96726-wl 

-2.6*096-03 

1.317*6-01 

339.03 

4.4*6-06 

UI.M 

t.***3£— ?1 

-2.3  6316—03 

1.46156-01 

356.99 

3.776-03 

llw.3* 

l. 7S2O6-01 

-7. *2*36-03 

1.2922E-01 

351.8* 

1.056-06 

112. 

9.*«*26-j2 

-2.*»116-03 

9.497*6-02 

351.51 

9.636-05 

u*.?r 

1. 76116-07 

-2.41746-01 

5.76646-02 

337.60 

1.726-03 

lll.ov 

1 .1*616-6  2 

-2.16**6-03 

1.56416-02 

351.30 

6.876-03 

111. or- 

— 2.91626-42 

-2.31566-4? 

2.92336-02 

114.3* 

2.456-06 

1*0.09 

-7.1*66;— 02 

-2.266:6-03 

7.35216-02 

181.77 

3.606-05 

142.10 

-l  .1*7*6— ;51 

-2.222’F-Ol 

1.1*766-01 

111.11 

1.2*6-05 

l**.:v 

-l. *1276-11 

-2.17776-01 

1.50296-01 

18C. 83 

7.666-05 

1*6.90 

-1.77*te-Cl 

-2.11466-03 

1.77  3  ?6-<>l 

110.69 

8. 716-06 

141.0? 

-1.9*Mt-4l 

-2.J9776-03 

1.9*326-01 

18b. 62 

2.766-05 

190.00 

-1.9409.;-'  1 

-2.95176—03 

1.99*06— 01 

18C. 39 

5.866-05 

192.00 

—1.9101*  — C 1 

-2.011*6-03 

1.910*1-01 

18b. 60 

7.70E-06 

1**.CC 

—1  ,6  9*96— 01 

-1.97116—1 ? 

1. *9316-01 

180.67 

4.876-05 

IS*. 00 

-1.1493F-U 

-1.93286-03 

1.3*9*6-01 

160.62 

9. 276-05 

lse.c? 

-1.63006-02 

-1.19496-03 

1.13206-02 

181.23 

3.976-03 

16C.90 

-3.1*106-02 

-i.ai73r-oi 

3.1*636-02 

183.38 

4.676-05 

1*2.?? 

3.3*5*6-42 

-1.62146-ul 

3. 36036-02 

356.69 

4.7*6-05 

1*4.0  . 

l.oif'r-r.i 

-1.7173J-03 

1.03116-01 

339.01 

9.906-09 

166.00 

1.76*26-01 

-1.73?  46-93 

1,76*36-01 

359.43 

3.966-05 

1*6. ?r 

2.*fo36-Ol 

-1.72036-03 

2 .  *8v36— >1 

359.60 

2.336-09 

17./. 4* 

l.ll*lf-vl 

-1.70(45-03 

3.13*16-01 

339.69 

5.976-05 

17>.'  ' 

1.7*42  -t  l 

-1. *76*6-03 

3.75*36-01 

339.7* 

7.7*6-05 

11* .v. 

*.?1’*6-01 

-1 .*60*'— 03 

4.253*6-01 

359.7* 

9.596-05 

176.  V 

*.*77*r-  l 

-i.O«7«E-o3 

4.62736-01 

35«.80 

9.00E-05 

17'. 

*.«‘97‘-..l 

-1.94f ":-43 

4.65976-01 

159.61 

.  9.026-05 

Fig.  8.  (Continued) 
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*  M  5  L  L  SO^ntt  »  A  0  I  A  l  VELOCITY 

FOOtlEH  TITLE  -  FOTNT-OOm*  5*Mt«6 

F606L6H  1.  -  GUTVA80  CONCENTMTEO  FOtCE  *T  NOtTH  POLE. 

Kt«  2.66 

FOEOUFttCV  iMt7.lt  HZ 


COUTITl'OEtCECYEESI 

VELOCITY 

A6S0LUTE  V4LUE 

FHASEIDE6) 

E6R06 

flac 

3.00 

4.35  476-6  3 

-2.fl4tf-01 

2.5144E-01 

270.99 

9.T2E— 05 

85 

2. J» 

4.J4Y1I-J1 

-2.3667E-01 

2.36915-01 

271.04 

9.405-05 

56 

♦  .nc 

4.29S6E-I  1 

-7.11136-Gl 

2.1119E-01 

271.17 

6.46E— 05 

34 

6.00 

4.2104r— ?3 

-1.73475-01 

1.7352E— 01 

271.39 

6.13E-05 

50 

P  *w 

4.C9325-C3 

-1.S3.6F-L  1 

1. 3312E-01 

271.76 

9.635-05 

39 

10. Qi' 

3.946CF-*!! 

-9.16*  65-02 

9. 1691c-02 

272.47 

3.86E— 05 

32 

12.00 

3.77136-03 

-J.2467F-J2 

5.2622E-42 

274.11 

4. 066—05 

42 

14.00 

!.!7!!KI 

-l.JM'f-M 

1.4T165-U2 

283.14 

6.056-03 
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The  velocity  calculation  also  deaonatrates  the  usefulness  of  the  data  parameter. 
NTERMS,  for  bounding  the  number  of  series  terms.  Although  the  smaller-magnltude 
results  for  colatitudes  38*  and  112*  are  not  fully  converged,  they  are  sufficiently 
converged  at  90  terms,  relative  to  neighboring  values  that  they  can  be  useful  with¬ 
out  further  refinement.  For  the  case  of  nonconvergence,  ERROR  could  have  been  set 
to  the  last  value  obtained  for  the  left  member  of  Eq.  41,  and  FLAG  equal  to  NTERMS. 
The  zero  values,  on  the  other  hand,  distinguish  more  clearly  such  cases  in  the 
column  listings. 

The  edit  sumaurlzing  polar  velocity  versus  ka  remains  from  an  earlier  version 
of  RADSPHERE,  before  the  general  capability  for  calculating  at  a  single  colatitude 
was  Introduced.  This  edit.  In  contrast  to  the  general  one,  treats  only  velocities 
at  colatitude  0*. 

To  plot  the  results  stored  internally  from  the  RADSPHERE  run,  the  additional 
data  in  Table  6  are  supplied  to  PLOTTER. 


Table  6.  PLOTTER  input  for  problem  of  point-driven  spherical  shell. 


Data  Item 

Value 

NKA 

1 

ORIG,  RMAX,  STEPSIZE 

0.,  13.,  2.6 

(for  surface  pressure) 

ORIG,  RMAX,  STEPSIZE 

-0.26,  0.26,  0.104 

(for  surface  velocity) 

- 

The  calculated  absolute  pressures  have  been  multiplied  by  the  shell  radius 
squared  for  comparison  with  the  published  results.  Figure  9  shows  the  echo  of 
identification  data  which  heads  the  data  file  being  plotted.  The  data  to  be  plotted 


on  this  file  can  also  be  echoed,  if  necessary,  by  using  the  flag  NDIAGN  in  the 


RADSPHERE  run. 


IIU  TO  M  PLOTTED  IS  FROH  THC  MOUM  RAOSPHERE 


PROQLEN  TIIU  -  POINT-DRIVEN  SRHCRI 

LOAD INC  CONDITION 

TO I NT  LOAD 

NACNITUOC  >  1.00 


RROOIEN  TYPES  THAT  HILL  Di  PLOTTED 


1 


QUANTITIES  AVAILADLE  FOR  PLOTTINC 


SHELL  SURFACE 
VELOCITT 


SHELL  SURFACE 
PRESSURE 


FAR  -  FIELD 
PRESSURE 


OATA  CALCULATED  AT  - 


INITIAL  COLATITUOt  •  0.00 
FINAL  COLATITUOS  >100.00 
OELTA  CO LATITUDE  ■  2.00 


OATA  IS  FOR  1  INTERVALS  IN  KA 


PLOTS  TO  IE  NAOS  ARE  POLAR  TTPE 


ENO  OF  DISSPLA  1.2  ~  *04*  VECTORS  SEMERATED  IN  2  PLOT  FR ARES. 

-ISSCO-  4116  SORRENTO  VALLEY  RLV0..SAN  D1E60  CALIF.  02121 

DISSPLA  IS  A  CONFIDENTIAL  PROPRIETARY  PRODUCT  OF  ISSCO  AND  ITS  USE 
IS  SURJECT  TO  A  NONOISSEHINATION  AND  NONDISCLOSURE  ASREENENT. 


Fig.  9.  Printer  output  from  PLOTTER  for  point-driven  spherical  shell. 


2 

The  RADSPHERE  results  for  a  |Surface  Pressure  |,  Fig.  10,  are  generally  in 

3 

good  agreement  with  Hayek's  results,  Fig.  11,  with  respect  to  profile.  To 
facilitate  comparison,  the  graph  in  Fig.  10  has  been  labeled  with  respect  to  the 
coordinate  system  used  for  computing. 

Figure  12  illustrates  the  method  (see  pages  19-20)  used  to  plot  shell  surface 


radial  velocities 


POINT-DRIVEN  SPHERE 
SHELL  SURFACE  RADIAL  VELOCITY 
POINT  DRIVE  -  MAG  -  LOO  (N.  POLE) 
FREQUENCY  -  MK7.322  Hz  ka  -2.063 


wm 

S«1 


Fig.  12.  Surface  radial  velocity  for  point-driven  shell, 
Q  -  0.85  and  h/a  “  0.03. 


SEC10R-DRIVEN  SPHERICAL  SHELL 

The  third  problem  is  from  a  set  of  calculations  used  in  the  validation  of 
NASHUA.  The  shell  is  driven  by  a  uniformly  distributed  load  centered  about  the 
pole,  0  ■  0°,  Fig.  2,  and  having  polar  angle  a  -  36°.  Given  the  load  magnitude, 

Fq  ■  20,  it  is  required  to  compute  the  field  pressure  at  radius  R  ■  100  for  ka  *  5. 
The  data  for  RADSPHERE  are  given  in  Table  7. 


K 


I 


Data  Item 

Value 

Title 

Sector-Driven  Sphere 

NFLAGF 

2 

Fo 

20. 

a 

36. 

a 

5. 

h 

0.15 

ps 

7669. 

V 

0.3 

E 

2.07E11 

n 

0. 

p fluid 

1000. 

c 

1524. 

NDIAGN 

1 

NPROB 

1 

NFLAG 

3 

NRADII 

1 

R 

100. 

NTERMS 

40 

e 

0.0001 

91 

0. 

02 

180. 

A0 

2. 

NINT 

1 

(ka1)1,  (ka^,  (Aka)  j 

5 . ,  5.,  1. 

PRO ALIN  Tint  -  SECTOR-DRIVEN  SPHERE 


DATA  FOR  SHELl 

RADIUS  *5.00 

THICKNESS  ■  ,13 

RAIS  DENSITY  >7.66401*03 

POISSON  RATIO  •  .30 

VOUNC  S  H0DULUS>2.0700E*11 
loss  FACTOR  *0.00 


DATA  FOR  EXCITATION 
SECTOR  LOAD 

POLAR  ANCLE  -  36.00  DECREES 
RA6NITU0E  ■  20.00 


DATA  FOR  FLUIC 

NASS  DENSITY  •1.0000006*03 
SPkEC  OF  SOUND  • 1. 52*000; *03 


OATA  FOR  F IS ID 


NO.  OF  RADII  ■  1 

RADIUS!  1)  >100.00 
RROtlEN  TYPES  TO  «E  CALCULATED 


1 


OUANTITIFS  TO  El  CONFUTED 

SHELL  SURFACE  SHELL  SURFACE  FAR  -  FIELD 

VELOCITY  PRESSURE  PRESSURE 

X 


DIAGNOSTICS 


DO  NOT  PRINT 


CALCULATIONS  TO  RE  HAOE  FOR  - 


INITIAL  COLATITUDE  ■  0.00 
FINAL  COLATITUDE  >160.00 
OELTA  C0LAT1TU0E  -  2.00 


Fig.  13.  RADSPHERE  output  for  sector-driven  spherical  shell 


CONV  Ef>G  t.NCc  C*IT£SIC,N  ■  1.0000E-04 


KA  INTERVAL  MiK8‘ A  i$  1 

INITIAL  A  A  ■  5. COOOL ♦ 00 
FINAL  KA  ■  5. OOOOfc ♦ 00 
OELTA  KA  -  1.00001*00 


Fig.  13.  (Continued) 
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»  I  1  l  0  MiSSUH  *  I  RADIUS  «  100.00 


’■Mien  mu  -  SICTOR-ORIVtN  SPM1R1 
PR06L1H  1.  -  0UTMR0  DISTRIBUTES  lORCt  AT  NORTH  ROLE. 
KA«  S.00 

■RtOUlNCV  •  262.5!  HZ 


Fig.  13.  (Continued) 
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6.89001-01 

356.26 
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12  Vf 
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12  £j 
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o 

6.66511-01 

-1.30571-02 
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12  £9 
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11  jCl 

It, 00 

6.3672L— 01 

6.6996E— 03 
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6.19351-01 
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12  Si 
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3. 9911c— 01 

3.28166-02 

6.00651-01 
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<j 

26L-95 

12 
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12  Sj 
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e 
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12 
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l 

011-05 

16  U| 

36.00 
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2.21311-01 

26.65 

3 

071-06 

16  K3 

3c. 00 

1,71121-01 
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1.9795E-01 
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5 

921-06 

16  S3 

6C.C0 

1.6690t-01 

6.03776-02 

1, 7516t-01 

36.17 
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971-05 

$1 
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12 
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12  vl 

66.00 
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16  ,>1 
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4*. 00 
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Data  for  PLOTTER  to  graph  these  results  are  in  Table  8. 


Table  8.  PLOTTER  input  for  problem  of  sector-driven  spherical  shell. 


Data  Item 

Value 

NKA 

1 

ORIG,  RMAX,  STEPS IZE 

0.,  0.60,  0.12 

The  problem  identification  echo  from  PLOTTER  is  shown  in  Fig.  14  and  the  plot 
in  Fig.  15. 

CAT*  TO  S;  PL'JTTrD  IS  FROM  THE  PROGRAM  PROS* HERE. 


PROSLFN  TITLE  -  SCCTOR-DRIVEN  SPNtRt 

LOADING  CONCITION 
SECTOR  LOAO 

POLAR  ANiU  >  36.00  CtCREtS 
P  AG  N I T  U  Dr.  ■  20.00 

PROBLEM  TYPES  THAT  MILL  6:  PL0TT10 

1 

QUANTITIES  AVAILABLE  EOS  PLOTTINL 

SHELL  SURFACE  SHELL  SURFACE  PAR  -  FIELD 

VELOCITY  PRESSURE  PRESSURE 

V 

NO.  OF  RADII  PCR  FI.Lt  CATA  Aft  •  l 
DATA  CALCULATED  AT  - 

INITIAL  COLATITUOE  ■  0.00 

‘INAL  COL  AT ITUDE  "11:0,00 
DELTA  C01ATITU0E  ■  2.00 

DATA  IS  FOR  1  INTEPVAlS  IN  KA 

PLOTS  Jj  Be  HADE  ARE  POLAR  TYPE 

ENC  OF  01SSPLA  8.2  —  3192  VECTORS  GENERATED  IN  1  PLOT  FRAMES* 

-ISSCO-  A1B6  SORRENTO  VALLEY  BLV0..SAN  DIEGO  CALIF.  92121 

01SSPIA  IS  A  CONFIDENTIAL  PROPRIETARY  PRODUCT  Of  ISSCO  ANO  ITS  USE 
IS  SUBJECT  TO  A  N ON D I SS EM  IN  AT  I OS  ANO  NONDISCLOSURE  AGREEMENT. 

Fig.  14.  Printer  output  from  PLOTTER  for  sector-driven  spherical  shell. 
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Fig.  15.  RADSPHERE  far-field  pressure  for  sector-driven  spherical  shell. 


Figure  16,  a  special  purpose  plot,  shows  superimposed  the  computed  results 


from  RADSPHERE  and  NASHUA  for  this  problem.  Agreement  is  quite  good. 
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FAR-FIELD  PRESSURE.  |p| 
AT  RADIUS  -  100 
ka-S 

— — —  RADSPHERE 
-  NASHUA 


10.  IS  1 0.30  To.46 


Fig.  16.  Comparison  of  RADSPHERE  and  NASHUA'S  results  for  the 
sector-driven  spherical  shell. 


RESONANCE  SEARCH 

In  the  course  of  validating  NASHUA*  for  ax i symmetric  vibration  of  submerged 
spherical  shells,  having  the  resonant  frequencies  of  both  structural  and  field 
response  was  essential,  so  that  those  ka  values  could  be  determined  at  which  the 
numerical  calculations  might  be  strongly  influenced  by  such  effects.  Our  fourth  and 
final  calculation  demonstrates  a  search  for  resonances  of  the  shell  used  in  the 
previous  calculation.  Data  to  run  the  search  were  obtained  by  modifying  the 


1 


data  (in  Table  7)  for  the  sector-driven  spherical  shell.  The  modified  data  for  this 
problem  are  shown  in  Table  9. 

Table  9.  RADSPHERE  input  for  resonance  search. 


Data  Item 

Value 

Title 

Resonance  Search 

NFLAGF 

• 

• 

Same 

• 

NPROB 

NFLAG 

NRADII 

5 

R 

Omit 

NTERMS 

60 

e 

61 

Same 

62 

0. 

A0 

1. 

NINT 

Same 

(ka^) l , 

(ka2)1,  (Aka) 1 

0.01,  10.,  0.005 

Figure  17  gives  the  data  input  echo  from  RADSPHERE,  followed  by  an  abbbreviated 
version  of  the  rather  lengthy  computed  results. 

The  appropriate  plot  for  this  type  of  calculation  is  y-log,  x-linear.  When 
the  data  for  RADSPHERE  specify  9}  -  02  the  program  assumes  that  response  versus  ka 
is  being  computed  and  sets  a  flag  for  log  plots  on  the  file  generated  for  PLOTTER. 

The  data  supplied  to  PLOTTER  are  given  in  Table  10. 


data  echo 


PRQQLER  TITLE  -  RESONANCE  SEARCH 


DATA  FOR  SHELL 

RAOIUS  -9*00 

THICKNESS  •  *19 

NASS  DENSITY  »T*66V0E«03 

POISSON  RATIO  •  .30 

VOUNC  $  H00ULUS-2.0T00EU1 

LOSS  FACTOR  >0.00 


OATA  FOR  EXCITATION 
SECTOR  LOAD 

POLAR  ANKLE  •  36*00  DECREES 
NASNITUOE  •  20*00 


OATA  FOR  FLUID 

NASS  DENSITY  -1*OOOOOOE«03 
SPEED  OF  SOU NO  •1*924000E*03 


PRQDLEN  TYPES  TO  RE  CALCULATED 
1 

FAR  -  FIELD 
PRESSURE 

0IA6NQSTICS 


QUANTITIES  TO  RE  CONPUTEO 


SHELL  SURFACE 
VELOCITY 


SHELL  SURFACE 
PRESSURE 


DO  NOT  PRINT 


CALCULATIONS  TO  RE  HADE  FOR  - 


INITIAL  COLATITUDE  •  0*00 
FINAL  COLATITUDE  •  0.00 
DELTA  COLATITUOE  •  1*00 


RAX*  NO*  OF  TERNS  FOR  A  SERIES  •  60 


CONVERQENCE  CRITERION  •  I.Q000E-04 


INTERVAL  NUNRER  IS  1 

INITIAL  KA  •  1*00001** 02 
FINAL  KA  •  1*0000E«01 
DELTA  KA  •  3.0000E-03 


.  RADSPHERE  output  for  resonance  search. 


Fig.  17 
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SHELL  SURFACE  RADIAL  VELOCITY 


COLATITUOE (DECREES I  ■  0.00 

DROOL EH  TITLE  -  RESONANCE  SEARCH 

RMOIEN  1 •  -  OUTNARD  OISTRIDUTED  FORCE  AT  NORTH  DOLE. 


RA 

FREQUENCY (HZ) 

VELOCITY 

AISOIUTE  VALUE 

FHAS2I02S) 

ERROR 

FLA6 

.010 

.409 

6.0000E-11 

2. 09602-04 

2.09602-04 

00.00 

7.942-09 

3 

•01S 

•  T2R 

1.3901E-10 

1.00372-04 

1.00372-04 

00.00 

3.742-09 

6 

.0X0 

.070 

2.4001E-10 

1.42722—04 

1.42722-04 

00.00 

6.662-09 

6 

.0X9 

1.213 

3.T40TS-10 

1. 14122—04 

1.14122-04 

00.00 

9.922-09 

7 

.0  JO 

1.499 

9.  3 00 Of  “10 

0.90412-09 

0.90412-09 

00.00 

7.962-09 

7 

.019 

1.000 

7.34761-10 

0.14702-09 

0.14702-09 

00.00 

0.402-09 

11 

.040 

1.040 

0.90932-10 

7. 12902-09 

7.12902-09 

00.00 

4.102-09 

12 

.049 

2.103 

1.2142E-00 

6.32022-09 

6.32022-09 

00.00 

9.302-09 

12 

•  090 

2.420 

1.40R7E— 00 

9.60212-09 

9.60212-09 

00.00 

6.992-09 

12 

•  099 

2.000 

1.R130E-00 

9.17042-09 

9. 1704E-09 

00.00 

7.032-09 

12 

.000 

2.011 

2.19712-00 

4.73932-09 

4.73932-09 

00.00 

0.492-09 

12 

.009 

3.193 

2.91002—00 

4.36262-09 

4. 36262-09 

00.00 

0.922-09 

16 

.070 

3.300 

2.0  3442-00 

4.04612-09 

4.04612-09 

00.00 

0.002-09 

16 

.079 

3. 030 

3. 36792-00 

3.77172-09 

3.77172-09 

00.00 

3.93E-09 

17 

•000 

3.001 

3.R302E— 00 

3.93102-09 

3.93102-09 

00.00 

4,402-09 

17 

•  009 

4.123 

4.32242-00 

3.31042-09 

3.31042-09 

60.00 

9,062-09 

17 

.000 

4.300 

4.R4402-00 

3.12022-09 

3.12022-09 

00.90 

9.602-09 

17 

•  009 

4.00E 

9.3090E-00 

2.09062-09 

2.09062-09 

09.90 

6.342-09 

17 

.100 

4.091 

9.07942-00 

2.00672-09 

2.00672-09 

00.00 

7,042-09 

17 

•  109 

9.004 

6.90402—00 

2.66012-09 

2.66012-09 

00.90 

7.702-09 

17 

.110 

9.336 

7.22372-00 

2.94102-09 

2.94102-09 

00.00 

0.992-09 

17 

.119 

9.9T0 

7.001 92-00 

2.42642-09 

2.42642-09 

09.90 

9.372-09 

17 

.120 

9.021 

0.9EE2E-00 

2.32271-09 

2.32272-09 

09.90 

9.602-09 

21 

.129 

6.064 

0.31302-00 

2.22402-09 

2.22402-09 

09.90 

3.992-09 

22 

.140 

6.406 

i.nnA6F-nA 

?.1446F>04 

2.I346F-04 

09.97 

3.642-04 

?? 

Fig.  17.  (Continued) 
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.119 

*>94* 

1.0*926-0* 

2.09076-09 

2.05076-05 

•9.97 

4.166-05 

•  140 

6.7*1 

1.16*91-0* 

1.47276-09 

1.47276-05 

•9.47 

4.496-05 

.149 

7.094 

1.29046-0* 

1. *9996-05 

1. *9996-05 

•4.96 

4. *16-05 

.190 

7.277 

1.99726-0* 

1. *31*6-05 

1. *11*6-05 

•9.96 

5.1*6-09 

.199 

7.91* 

1.42696-0* 

1.76796-05 

1.76796-05 

■9.95 

5.546-05 

.140 

7.7*2 

1.91*96-0* 

1.707*6-05 

1.707(6-05 

•9.95 

5.926-05 

•  1*9 

0.004 

1.614*6-0* 

1.69136-09 

1.69136-05 

•9.94 

6.326-05 

.170 
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1.71296—00 
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1*59746-05 
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.179 
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1.54746-05 
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7.156-05 

.100 

1.792 
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1.49966-05 

1.49966-05 

•9.93 

7.546-05 

.109 

6.974 

2.021*6-0* 

1.49426-05 

1.45426-05 

•9.92 

*•056-05 

.1*0 

9.217 

2.19906-0* 

1.41116-05 

1.41116-05 

•9.91 

(.526-05 

.1*9 

9.460 

2.24496-0* 

1.17016-05 

1.17016-05 

•9.91 

9.006-05 

.200 

4.702 

2.99996-0* 

1.31106-09 

1.33106-05 

•9.90 

9.506-05 

•  209 

«.*49 

2.476*6-0* 

1.292  96-09 

1.29216-05 

•4, *9 

1.936-05 

.210 

10.167 

2.59676-0* 

1.25666-05 

1.25666-05 

•  9.99 

4.146-05 

.219 

10.490 

2.  71446-08 

1.22256-05 

1.22256-05 

*9.97 

4.356-05 

•  220 

10.672 

2.04476-08 

1.1*976-09 

1.1*976-05 

(9. *6 

4.5*6-05 

.229 

10.419 

2.97266-0* 

1.15(46-05 

1.15*46-05 

•9. *5 

4.916-05 

.290 

11.197 

9.10916-0* 

1.12*26-05 

1.12*26-05 

(9.(4 

5.056-05 

.299 

11.400 

9.29616-0* 

1.09496-05 

1.09916-05 

•4. *3 

5.246-05 

.240 

11.649 

9.97206-0* 

1.07146-09 

1.07146-05 

*9.92 

5.556-05 

.249 

11.6*9 

1.51096-0* 
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Table  10.  PLOTTER  input  for  resonance  search. 


Data  Item 

Value 

NKA 

1999 

NCYCLE,  YMIN 

7,  l.E-7 

XSTEP 

1. 

NXT1CK,  NYTICK 

-2,  -9 

Figure  18  gives  the  echo  from  PLOTTER  identifying  the  data,  and  Fig.  19  shows 
the  response  curve  of  radial  velocity  at  0  *  0°  versus  ka,  whose  peaks  indicate 
the  resonances. 

For  this  calculation  the  total  ka  computed  is  just  two  under  the  maximum 
currently  allowed  by  the  program.  In  the  event  that  more  resolution  is  needed,  the 
primary  Interval  can  be  subdivided  and  the  same  number  of  points  used  over  each 
•ublnterval.  Computation  over  the  entire  set  of  subintervals  is  done  in  a  single 
run  of  RADSPHERE. 

The  1999  points  in  the  velocity  response  curve  were  calculated  in  9.462  CP 


seconds  on  the  Center's  CDC  176  computer. 


DATA  TO  BE  PLOTTED  IS  FROM  THE  PROGRAM  RADSPHERE. 

PROBLEM  TITLE  -  RESONANCE  SEARCH 

LOADING  CONDITION 
SECTOR  LOAO 

POLAR  ANGLE  •  36.00  OEGREES 
MAGNITUDE  ■  20.00 

PROBLEM  TYPES  THAT  MILL  BE  PLOTTED 

1 

QUANTITIES  AVAILABLE  FOR  PL0TTIN6 

shell  surface  shell  surface 

VELOCITY  PRESSURE 

X 

DATA  CALCULATED  AT  - 

INITIAL  COLATITUOE  •  0.00 

FINAL  COLATITUOE  •  0.00 

DELTA  COLATITUOE  »  1.00 

DATA  IS  FOR  1  INTERVALS  IN  KA 

PLOTS  TO  BE  MADE  ARE  LOG-LINEAR  TYPE 

END  OF  OISSPLA  B.2  —  3672  VECTORS  GENERATED  IN  1  PLOT  FRAMES. 

-ISSCO-  41R6  SORRENTO  VALLEY  8LVD..SAN  OIEGO  CALIF.  92121 

OISSPLA  IS  A  CONFIDENTIAL  PROPRIETARY  PROOUCT  OF  ISSCO  AND  ITS  USE 
IS  SUBJECT  TO  A  NONDISSEMINATION  AND  NONDISCLOSURE  AGREEMENT. 

Fig.  18.  Printer  output  from  PLOTTER  for  resonance  search. 
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1 .  DTNSRDC  reports,  a  formal  aortas,  contain  information  of  permanent  technical 
value.  They  carry  a  consecutive  numerical  identification  regardless  of  their  classification 
or  the  originating  department. 

2.  Departmental  reports,  a  semiformal  series,  contain  information  of  a  preliminary, 
temporary,  or  proprietary  nature  or  of  limited  interest  or  significance.  They  carry  a 
departmental  alphanumerical  identification. 

3.  Technical  memoranda,  an  Informal  series,  contain  technical  documentation  of 
limited  use  and  interest.  They  are  primarily  working  papers  intended  for  internal  use. 
They  carry  an  identifying  number  which  indicates  their  type  and  the  numerical  code  of 
the  originating  department.  Any  distribution  outside  DTNSRDC  must  be  approved  by 
the  head  of  the  originating  department  on  a  case-by-case  basis. 
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